Subleading Corrections To Thrust Using Effective Field Theory 



Simon M. Freedman 1 

1 Department of Physics, University of Toronto, 60 St. George Street, Toronto, Ontario, Canada M5S 1A7 

(Dated: March 8, 2013) 

We calculate the subleading corrections to the thrust rate using Soft-Collinear Effective Theory to 
factorize the rate and match onto jet and soft operators that describe the degrees of freedom of the 
relevant scales. We work in the perturbative regime where all the scales are well above Aqcd. The 
thrust rate involves an incomplete sum over final states that is enforced by a measurement operator. 
Subleading corrections require matching onto not only the higher dimensional dijet operators, but 
also matching onto subleading measurement operators in the effective theory. We explicitly show 
how to factorize the 0(a s r) thrust rate into a hard function multiplied by the convolution of the 
vacuum expectation value of jet and soft operators. Our approach can be generalized to other jet 
shapes and rates. 



I. INTRODUCTION 

Jet shapes are examples of observables with multiple 
scales, which can give rise to logarithmic enhancements in 
the fixed order rate that ruin the perturbative expansion 
in the strong coupling constant. Effective Field Theo- 
ries (EFTs) separate the scales by expanding in their 
small ratio and matching onto operators that describe 
the degrees of freedom at each scale. The logarithmic 
enhancements are then summed by using the renormal- 
ization group to run between the scales restoring pertur- 
bative control of the rate. The effects of the subleading 
corrections from the small ratio of scales are systemat- 
ically calculated in the EFT by matching onto higher 
dimensional operators. 

The appropriate theory for describing jet shapes with 
narrow jets is Soft-Collinear Effective Theory (SCET) 
P~H7]. Deriving factorization theorems that separate the 
scales is straightforward at leading order (LO) in SCET 
due to the explicit decoupling of the collinear and soft de- 
grees of freedom. Subleading corrections to SCET have 
previously been used to study subleading corrections to B 
decays [Ej; however, jet shapes are complicated by kine- 
matic cuts placed on the phase space. The goal of this 
paper is to demonstrate how to factorize the subleading 
corrections to jet shapes using SCET by considering the 
example of the thrust observable. 

Thrust [9] is defined by 
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T = — 22 mil1 ( E * + **' E i ~ Pi) 
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where Q is the initial energy, X is the final state, and the 
thrust axis t is the unit vector that maximizes ^2 ieX K ' 
Pi | . When t< 1 the final state is a pair of back-to-back 
jets with small invariant mass. Thrust is a convenient 
observable to illustrate how to calculate subleading jet 
shapes in SCET due to its simple phase space. We will 
concern ourselves only with e + e~ "/* — > X in order 
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to reduce the contribution from initial state radiation 
and restrict ourselves to vector currents. It is possible to 
generalize to weak process but for simplicity, we will not 
consider axial currents in this paper. 

The LO in t thrust rate was calculated using SCET in 
[TCTI ITT] . The rate was written in the factorized form 

H x (0| J|0) <g> (0| J|0) ® (0\S\0) + 0(t), (2) 

where the convolution is defined in Sec. IIHI The LO rate 
was factorized by matching onto SCET and expanding 
the final state phase space imposed by ([T]) in the SCET 
power counting. The jet operators, J and J, describe the 
physics at the intermediate scale \frQ, the soft operator, 
S, describes the physics at the soft scale tQ, and the hard 
function, H, describes the physics above the hard scale 
Q. Factorization allows the large bar's to be summed by 
separately renormalizing the jet and soft operators. This 
was done in 

We are interested in extending the results of [TOJ [TTJ to 
include the 0(t) corrections to the rate. We restrict our 
calculation to the perturbative regime when the soft scale 
is well above Aqcd ■ We will ignore the effects of hadron 
masses, which have been discussed in jJ5]. We will show 
how the 0(t) rate can be factorized as in ^ using SCET 
with the appropriate subleading jet and soft operators. 
The se subleading operators are generalizations of the LO 
operators, and properly separate the scales while having 
consistent power counting. We leave renormalizing these 
operators to a future work. 

Understanding how to incorporate subleading phase 
space effects is important for writing a factorization the- 
orem beyond LO. We will show how subleading phase 
space is accounted for in the effective theory by 1) consis- 
tent expansion of the cuts using the SCET power count- 
ing, and 2) insertions of subleading operators that ac- 
count for the QCD momentum conservation expansion. 
Both these effects must be accounted for in order to cal- 
culate the 0(t) corrections to the rate and reproduce the 
perturbative QCD result at 0(a s r) in [T5] . 

The rest of the paper is organized as follows: in Sec. [TTJ 
we review the description of SCET in [7] and explain our 
reasoning for using this formulation over the formulation 
in [THE]. In Sec. Ill we review the LO calculation and 
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introduce the notation used in the rest of the paper. We 
calculate the 0(a s r) rate in Sec. IV and demonstrate how 
to write it in a factorized form. We conclude in Sec. |V] 
The full list of operators is reserved for the appendices. 
For simplicity we use the notation Lh = ln(/i 2 /Q 2 ), 



L., = \n(^/(Q^Y), L t 



a s C F /{2ir) where Cp 



(N 2 C 



1u(ij, 2 /(Qt) 2 ), and a s = 
■1)/(2N C ) fovN c colours. 





II. REVIEW OF SCET 



FIG. 1: QCD vertex diagrams required for dijet 
production at 0(a s ). 



Soft-Collincar Effective Theory describes the interac- 
tions between highly boosted and low energy degrees of 
freedom. Three types of fields are required for calculat- 
ing thrust: n-collinear, n-collinear, and soft, which have 
characteristic momentum scaling in terms of the SCET 
expansion parameter A< 1 



K~Q(A 2 ,1,A), 



p£~Q(l,A 2 ,A), 



k^QX 2 (3) 



respectively. The two collinear sectors describe fields 
that are highly energetic and moving in opposite direc- 
tions, while the soft sector describes low energy fields. 
The SCET Lagrangian and operators are derived by ex- 
panding the interactions between the different fields in A. 
Momenta denoted by k will be 0(A 2 )Q unless otherwise 
specified. 

We will use the SCET formulation of [7] . In this formu- 
lation, soft and collinear fields are described by QCD in 
the absence of an external current. Therefore, the SCET 
Lagrangian is 



r — r n A- r n a- r s 

SCET — *~QCD ' QCD ~ QCD 5 



(4) 



which has no subleading contributions and each £™ CD 
will describe fields from the m sector only. The interac- 
tions between the sectors are contained in the external 
currents. The QCD vector current 

J 2 "{x) = $(x)-fil>(x) (5) 

is matched onto the SCET dijet operators 

c£ ) 0<£ ) > 1 + i y C^O^ + ... (6) 



i>i 



where the superscripts denote the suppression in A and 
the ellipses represent higher dimensional operators. The 
operators and their tree-level matching coefficients 
are found by expanding the diagrams in Fig. [I] in A for 
n-collinear, n-collinear, and soft fields. 

At LO, the (anti-)quark must be (n-)n-collinear and 
the gluon can be cither collinear or soft. The LO dijet 



1 The momentum p M = (p + , p , p± ) is denned in light-cone coor- 



dinates by jv 
Pi = P M " P 4 



p ■ n 



E — p ■ n, p 
t- where n ■ n ■ 



p ■ n = E + p ■ n, and 



operator is [7] 



(7) 



with one-loop matching coefficient [M] 
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and MS counterterm 



Z { °\a) =l + a s 



0{a 2 s 



3 
2~c 



(8) 
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-Q 2 



+ 0{a 2 s ). 
(9) 



The interactions between the different sectors are repro- 
duced by Wilson lines defined in a representation R 



F4 R >(x,y) = Pexp 
YW(x,y) = Pexp 



-w 



dsn- A^(x + ns)T^j 

dsn- A a s (x + ns)T^j 
(10) 



with similar definitions for and Y£ . The projec- 

tors P n = {i/ifi) 1 4 and P„ — {^aji) /4 are required from 
the expansion of the Dirac equation. The positions are 
defined as 

x n = (0,2! • n,Sx) x ^ = (0, oo,xj_) x^ s = (0, oo, 0) 
x fL = (x ■ n, 0, x±) x°° — (oo, 0, x±) x^ s — (oo, 0, 0) 

(11) 

and are necessary to conserve the appropriate compo- 
nents of momentum that respect Each square 
bracket in Q is a separately gauge invariant piece, which 
means the sectors explicitly decouple from one another. 
The physical interpretation of the operator is given in [7] . 

The next-to-leading order (NLO) dijet operators and 
matching coefficients were found in 7] and are repro- 
duced in Appendix]^] The operators are generalizations 
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of the LO operator ([7]) with appropriate derivative inser- 
tions. For 0(r) corrections to thrust, we will also need 
the N 2 LO operators, which are found by following the 
approach of [TJ. These operators and their matching co- 
efficients are also shown in Appendix [XJ Both the NLO 
and N 2 LO dijet operators explicitly decouple the sectors 
in the same way as the LO operator. 

A more widely used formulation of SCET [1-5J sep- 
arates the collinear momentum into p = p + k, where 
k ~ X 2 Q and p ~ Q, XQ are the residual and label mo- 
mentum respectively. The large label momentum is re- 
moved from all interactions leading to the Lagrangian be- 
ing a non-local expansion of two-component spinors with 
soft and collinear fields explicitly interacting at NLO. 
Although, the formalisms of [IH5] and [7J are equivalent, 
they approach momentum conservation differently. This 
is important when cuts are placed on phase space such as 
for the thrust rate. In the approach of [7j, momentum is 

not conserved and the subleading 0% operators in Ap- 
pendix [A] account for the expansion of QCD momentum 
conservation using SCET momentum power counting. In 
the approach of [IH5], label and residual momentum are 
separately conserved meaning momentum is exactly con- 
served. The equivalent action of the 0^ S ^ operator comes 
from the subleading kinetic interaction |15) 

4e L=o = J2^A^^t ± ^(x) + h.c (12) 

p,p' 

where t; n ,p is a two-component n-collinear spinor with 
label momentum p, and T" 1 is an operator that pulls 
down label momentum. Insertions of this term into time- 
ordered products is equivalent to expanding the on-shell 
condition (p + k) 2 = in A. These terms are not neces- 
sary for inclusive phase space calculations such as sub- 
leading B decays [5]. However, they will be important 
when cuts are placed on phase space such as in the sub- 
leading thrust rate and reproduce the action of the 2 U ^ 
operator in time-ordered products. 

In this paper we choose to use the formalism of [7J . The 
Lagrangian of [7J has simpler Feynman rules and only 

one insertion of 2 1S ^ is required instead of an insertion 
of £tl' for each collinear field. The explicit decoupling of 
sectors at the operator level also makes it easier to derive 
a subleading factorization theorem. 



III. LEADING ORDER CALCULATION 

The LO thrust distribution was calculated using the 
approach of [TH5] in [TT] and was written in the factorized 
form ([2|. In this section we review the calculation using 
the formalism of [7J , which gives an equivalent form of the 
answer. In the next section we generalize this description 
to calculate the 0(t) rate. 



The thrust rate is the cumulate of the distribution 

= J d^xe-^iOUi^ (x)M QCD (r) (0) |0) 

where — (Q/2)(n M + n^) is the momentum of the 
incoming photon and oo is the Born cross-section. The 
measurement operator, M QC d{t) PHIHE], acts 011 states 
\X) 

M qcd (t)\X) = M qcd (t, {px})\X) (14) 

to project only those final states that give a thrust value 
r. When taking the cuts of diagrams, the function 
■M- Q cv{t, {px}) generates the appropriate phase space by 
restricting the momentum of the particles {px}- 

To factorize the rate, we first match the QCD currents 
and measurement operators onto SCET dijet and mea- 
surement operators. The matching of the QCD currents 
onto SCET dijet operators was discussed in Sec. [TTJ The 
QCD measurement operators are matched onto SCET 
measurement operators in a similar manner 

AWr) A?°»(t) + M«(t) + M {2 \t) + 0(A 3 ), 

(15) 

where the superscripts refer to the suppression in A. The 
SCET measurement operators are found by expanding 
the thrust constraints implemented by A^qcd using the 
SCET momentum scaling ([3]). 

Thrust is measured with respect to the thrust axis, 
t, defined below |l]). The definition of the thrust axis 
in SCET has an expansion in A and is written as t = 
_|_ 0(A 2 ). The SCET momentum power counting 
defines the LO thrust axis = —n [TU], where we 
have chosen the — n axis to be exactly along the total 
n-collinear momentum (i.e. jpki = 0). The overall sign 
of t is unimportant as seen in (Jit) . The sectors decouple 
in the LO measurement operator [TUl E] 

( T ) = (r„) ® (r fi ) ® (r s ) (16) 

because the thrust axis is independent of any individual 
particle. The convolution above is defined as 

/i(n) 8/2(75)8/3(7*)= (17) 

J dTldT 2 dT 3 6(T -T1-T2- T 3 )fl(Tl)f 2 {T 2 )f3{T 3 ). 

Using the LO definition of the thrust axis, the action of 
the measurement operators are 

^<,, W ,-(^p( T -l£!W!) 

Mf\r, {k}) = 5 {r-^(n- kf t + n ■ k W t )^j (18) 
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where the sums are only over the momentum in each 
sector. We have defined 



0(A 4 ) (19) 



as the total soft momentum in the ±i hemisphere. The 
LO definition is 

k { T = £ K6(±h ■ m) = ]T K9( T h ■ ft), (20) 



where in both (191 and ( 20 ) the sum is over all soft parti- 
cles. The d — 4 — 2e dependent prefactors in the collinear 
sectors come from choosing the collinear fields to be in 
the n M and n M directions. At LO these prefactor have no 
affect, but are important for the 0(t) corrections. 

By matching the QCD operators in ( 13 1 onto the 



SCET operators, the LO thrust rate is written as 



R(t) = \C, 



_ ir»(o)i2 



d d 3 



(21) 
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(x)M {0) (t)0? V(0) |0) +0( 



(°), 



where the O(r) corrections will be calculated by the sub- 
leading in A operators. The rate is factorized by matching 
above the operator product onto jet and soft operators 



(*)M<°>(t)O<J(0) 



(22) 



with matching coefficient C^°K As usual, the superscripts 
on the jet and soft operators refer to their suppression in 
A. The rate can then be written in the desired factorized 
form 

R(T) = H^^)(0\J^(^r n )\0) (23) 
®(0|j(°)(/x,rn)|0) ® (0|5( Q )(M,r s )|0) + 0(r) 

where the hard function is the product of the matching 
coefficients, H^(ji) = |Cf 5 \n)\ 2 C^ (//). 

The explicit decoupling of n-collinear, n-collinear, and 
soft degrees of freedom in the dijet and measurement op- 
erators makes finding the appropriate jet and soft opera- 
tors straightforward. Using the Fierz identity to separate 
the spin and colour indices, we find the operators are 



/'■ D - ^-Tr | dz+d d -Ve-^„(0)W( 3 )^ 

^Tr | d*"^ Ve^*" (x n , O | A^ 0) (r)^i 3) (a^ , 0)^(0) 

-^Tr ri 3) , o) ri 3 ) (o, )M<°> (t)y^ (x~- , o)yi 3) (0, ) , 

JVC 



(24) 



where the trace is over spins and colour. These operators give the same Feynman rules as those found in |llj . The 
matching coefficient in ( 22 ) is most easily found by comparing the vacuum expectation value of both sides. The real 
emission contributions to the one-loop vacuum expectation value of and are pictured in Fig. [2] They are 
calculated by cutting the diagrams along the vertical dashed lines and applying the measurement operator to the 
fields passing through this cut [TT]. The virtual diagrams are scaleless and thus zero in MS. The vacuum expectation 
value of the jet and soft operators are [11] 



<0|J( V,t)|0)=<5(t) 



1 + a s 



-L 



H 



2e 



+ L 



H 



-L 



7 7T- 



H 



2L 



<0|j( )(M ! r)|0> = (0|j(°)( / x,r)|0> 



(0\SW(h,t)\0)=6(t) 1 - n. s 



(25) 



2 2 

72 ~ 7 



Tl- 



4L.c 



(r) 



where we have included the zero-bin procedure |17j . 
which accounts for the double counting between the 
collinear and soft operators. The matching coefficient 
C(°)(/z) = 1 [11] meaning the hard function is 

JT<°>(M) = l + a s (l 2 h - 3L H - 8 + ^ \ + 0(a 2 s ). 

(26) 

As expected, the jet and soft operators separate the y/rQ 
and tQ scales. The hard function describes the physics 



above the cut-off Q, of the effective theory. 

Matrix elements of the jet and soft operators also do 
not need any further expansion in t. This is required to 
ensure that operators of different orders do not mix. We 
note that this is different than the results in 18], which 
considered the exclusive JADE two-jet rate at LO. This 
rate has the same 0(a s ) phase space as the thrust rate 
in QCD. In [18], the phase space was not consistently 
expanded in A and an expansion in r was required af- 
ter the phase space integration. A subleading zero-bin 
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,® 



(e) (0|S<°)|0) 



® 



/x 



(f) (0|S(°)| 



® ® : ,® 



(g) (0[St°)[0> 



Xj X 

Y 

(h) (0|5W|0) 



FIG. 2: One- loop diagrams of (25 1. Solid lines represent fermions, dashed lines represent Wilson lines with the 
colour flowing in the direction of the arrows, and the dots represent Lagrangian insertions. The type of soft Wilson 



lines are labelled in Fig. 2e The cut is distinguished by the bold vertical dashed line. The contributions to look 



identical to the contributions of after a rotation of 180° 



procedure is necessary to ensure that the LO operators 
do not contribute to subleading corrections. In this pa- 
per and in , the measurement operator is consistently 
expanded in A so matrix elements of the jet and soft op- 
erators automatically have consistent power counting. 

The thrust rate at O(a s T ) is calculated by substitut- 
ing the hard function and the vacuum expectation values 
of the jet and soft operators into (23). The rate is found 
to be 



R(t) = 1 



-21n 2 r- 31nr 



- 1 



0(a 2 s 



(27) 

which reproduces the rate found in perturbative QCD at 
this order By separately renormalizing the jet and 
soft operators the large lnr's were summed in [11]. In 
the next section we follow the same procedure to find the 
0(a s T) correction to the rate and write it in a factorized 
form analogous to (23). 



IV. 



NEXT-TO-LEADING ORDER 
CALCULATION 



The results are extended to include the O(r) correc- 
tions by systematically matching the QCD current and 
measurement operator in ( 13 ) onto subleading SCET di- 



jet and measurement operators. The 0(t) thrust rate in 
SCET is written as 



(28) 



i+j + k<2 



x(0\O { 2 ^(x)M^(t)O { £(0)\0) + O(t 2 ) 



From the LO calculation A ~ y/r, so we need the 0(A 2 ) 
SCET operators as illustrated by the constraints on the 
sum. One can explicitly check that the O(A) corrections, 
which would give 0(y /f r) corrections, vanish. 



As in the previous section, we want to write ( 28 1 in 



a factorized form by matching onto subleading jet and 
soft operators. These operators will be generalizations of 
the LO operators and their matrix elements must have 
a consistent power counting in r as in the LO case. We 
will only do the tree-level matching, which is all that is 
necessary to calculate the 0(a s r) rate. The subleading 
dijet operators are written in Appendix [A] and explicitly 
decouple the sectors. We must also find the sublead- 



ing measurement operators, A4^' 2 \t) in (15). We will 
show in the next section that the action of the subleading 
measurement operators also decouples the sectors, anal- 
ogously to ( 16 ). The explicit decoupling of the sectors in 



the dijet and measurement operators makes writing the 
rate in the desired factorized form in Sec. |IVB| straight- 
forward. 



A. Measurement Operator 

The action of the subleading measurement operators 
are found by first expanding the definition of the thrust 
axis and then finding this expansion's effect on the mea- 
surement of thrust. The thrust axis is defined as the 
unit vector that maximizes the sum below 0. The sum 
is maximized when the thrust axis is in the direction of 
the hemisphere with the largest three-momentum, p +t 
[LO] . Therefore, the thrust axis is written in SCET as 
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the expansion 



r = £t* = r(°) + + + (A 6 ), 
\p+t\ 



(29) 



where the superscripts refer to the suppression in A. 

In order to find p+t, we first note that the n-collinear 
and n-collinear particles are always in opposite hemi- 
spheres. SCET momentum power counting enforces this 
at LO and the zero-bin procedure will enforce this at 
all orders in A. Therefore, the total momentum of the 
hemisphere will be 



p+t =Pn + k +t , 



(30) 



where p£ is the total n-collinear momentum and k^_ t is 
the total soft momentum in the +t hemisphere as defined 



in (19). The expansion of (19 1 and (29) allows us to 



iteratively solve for both f W and k^Y . The subleading 



corrections to the thrust axes are 



£(2) _ Zfi +t± 



pi 



(31) 



(pi 



J (4) _ 



:(o) 



+ ^2 



(pi) 2 



fc (0) 



2k 



(2) 
+t_L 



Pn 



where k~^ is the total n-collinear momentum in the n M 
direction. The subleading correction to the total soft 
momentum in each hemisphere 



•„(2)j" 

"+t 



,(2)M 



1 



■ fc«x) S(-n ■ ki) 



(32) 



is found by inserting the t ^ into ( 19 ). The first equality 



above is because the sum of the soft momentum in the 
two hemispheres must be 0(X 2 ). 

We note we could have instead used —p+t — Pn+k-t in 



( 29 1 , where p n is the total n-collinear momentum and k_ t 
is the total soft momentum in the —t hemisphere. How- 
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ever, the defi nition in ( 30 ) is simpler due to our choice of 
n in Sec. 



such that Pnj_ = 0. The apparent asymmc- 



space will be accounted for by the and O^ 5 ^ op- 

erators in Appendix |a| The choice of pn± = means the 
Feynman rules of these operators involve d/dp n ± , s only. 

The subleading measurement operators are found by 
substituting the corrections to the thrust axis into 
We first consider the contribution from an n-collinear 
particle with momentum p^. As discussed in the second 
paragraph of this section, Ei+pi-t is always the minimum 
for each n-collinear particle. Therefore, the n-collinear 
sector contributes 



+p l -t (o) ) + ^[t^ +t 



Q 



'(2) 



0(A 5 ) 

to the thrust, where the sum is only over n-collinear par- 
ticles. The first term reproduces the action of the LO 
measurement operator in (18). The term giv es the 
NLO correction and the first term of t*- 4 - 1 in (31) gives 
the N 2 LO correction. This power counting is due to 
p±i ~ 0(X)Q for n-collinear particles. 

The contribution to thrust from the n-collinear parti- 
cles is found in a similar way. Here, the Ei—pi-t is always 
the minimum so the n-collinear sector contributes 



s2> 



Q 



(34) 



where the sum is only over n-collinear particles. The first 
term is the LO contribution and the second term is the 
N 2 LO contribution. There is no i"^ term because we 
have set Pn± = by our choice of n. 

Unlike the collinear particles, soft particles can be in ei- 
ther hemisphere. The minimum of Ei±t-ki is determined 
by which hemisphere the soft particle is in. Therefore, 
the soft sector contributes 



1 

Q 



■ k$ +h-k 



(0) 



1 

Q 



„(2) 



(35) 

(m-m))+o(x% 



to the total thrust. The first line is the LO contribution 
in (18) and the remaining terms are all N 2 LO. 



The action of the subleading measurement operators 
is found by Taylor expanding the contribution to thrust 
in A. We incorporate the NLO and N 2 LO corrections 
from the collinear sectors in ( 33 1 and ( 34 ) by writing the 



try in the labelling of n^ 1 and n M in the resulting phase action of the subleading measurement operators as 



M [ln \T,{p n ,pM) = 

M {2n «\T,{Pn,Pn,k s }) = 



pZ±A 
Q dr, 



Pn±Pn± 92 

Q 2 dr 2 



M^(T n ,{p n }) ) ® ( ^- M g\rn,{Pn}) 



'91.(0)0, (0)(8 



7 



M i2nb) (T,{ Pn , Pn ,k s }) 



M( 2 ^)(r, {p n ,p n , k s }) = M^(r n , { Pn }) 



Mf{T ni { Pfl }) 



i^LM«°)(r a , {*.})) (36) 



Q 9t r 



_\^(,2n a ) comes f r0 m the sec ond term in the expansion of 
the NLO correction in (33). The N 2 LO corrections from 



the soft sector in ( 35 ) are incorporated by the subleading 



measurement operators 

M ( - 2si) (T,{p n ,p n ,k s })= (37) 

MW(T n ,{ Pn })® ®M^(T s ,{k s }) 

where 



M 2si) MM) 

M^(T,{k s }) 

M^(r,{k s }) 
M^(r,{k s }) 



2k 



(0) 



:(o) 



+t_L n '+t_L 



^2/c 



Q 

(o) r(o) 



dr. 



+t-L 



-t_L 



Q 

ptn ■ fe (2) 



M 0) (^{M) 
■^i 0) (^,{fcJ) 



+/ 



Q 2 



Or, 



Q 2 3r s 



■Mi°)(r s ,{fc s }) (38) 
M 0) (MU)- 



The p^" in the action of the soft measureme nt o perators 
j^[(2s ± t) are introduced to cancel the p% in (32 1. 



The actions of the measurement operators (36) and 
(37) define the NLO and N 2 LO measurement operators 



M( 1,2 )(t). As the brackets suggest, the sectors explicitly 
decouple in the action of the subleading measurement op- 
erators. This is due to the corrections to thrust depend- 
ing only on the total momentum of each collinear sector 
and not any individual particle. While it is possible to 
formally write the subleading measurement operators us- 
ing the energy-flow operator [TUl E] and not just their 
actions in momentum space, we see no reason to do so: 



only their actions in momentum space are necessary in 
calculations. 

We note that the measurement operators in this sec- 
tion are found using the formalism of [7] and would be 
different if we used the SCET formalism of [T]-[5]. It was 
suggested in [10] the subleading measurement operators 
could be found using the subleading terms in the SCET 
Lagrangian of [TH5]. While we do not explicitly check 
this, we note that the breaking of the explicit decoupling 
of soft and collinear fields in the subleading Lagrangian 
would complicate factorization. 



B. Factorization 

The explicit decoupling of the sectors in the subleading 
dijet and measurement operators makes it straightfor- 
ward to factorize the subleading corrections to the rate. 
In order to factorize the rate, each operator product in 
(28) is matched onto the appropriate jet and soft opera- 



tors 

d d xOf\x)M (k) (T)0 ( i\$) (39) 

a(l,m,n) 

+TC (i,j,k) ji0){Tn) 3 J( o )(Tfi) g, 5 (o )(Ts); 

with matching coefficients C^' k \ The last line is made 
N 2 LO by the explicit r in front. The integer a > 1 labels 
the combination of jet and soft operators and we require 
i + j + k = 2 = l + m + n. The 0(t) rate can then be 
written in the factorized form 



(r) - Y, M <°l J(l) T ») 1°) ® <°l J(m) (/*> T ») 1°) ® (A*, t.) |0> 



+T J2H^- k \fi)(0\J^(^T n )\0) ® <0|J< )(a*,7«)|0) ® (0|5(°)( M ,r s )|0) 



(40) 



where the hard functions arc defined as 



H 



a,0 



(m) 



•<(-t,J,fe)/ 



(41) 



This generalizes the LO factorization (23) to incorporate 
the N 2 LO corrections and is the main result of our paper. 



from (28), the factorized form will allow the jet and soft 



operators to be renormalized separately. 

Below we will use a few examples to demonstrate how 
the jet and soft operators in (39) are found. The full list 
of operators and their matching coefficients are found 



While it is possible to calculate the 0(a s r) rate directly in Appendix [B] For the sake of brevity, we only write 



those operators that contribute to the 0(a s r) rate. We 
omit the phase space integrals when calculating matrix 
elements of the operators to avoid overbearing the reader. 

A) i — la n , j — lb n , k — : The operators are 
found in ((Al. The left-hand side of ([391) is 



.,— iQ-x 



1 

Q 



$n{x n )W£\x n ,X™) 

x yi 3) (^,o)ri 3 '(o,^) 

^f'W ® A^ 0) (r B ) ® (r.) 

xi [^(0^(0)^(0,0' 

x '^(c.ojyf (o,a£")l M 8) OC*, 0)^(0) 



(42) 



where = f 7^7 M P n and = ^l^P n and 7^ is 
defined in Appendix |B| The 1/Q's come from the defini- 
tion of the operator expansion (JgJ) and the negative sign 
comes from taking the hcrmitian conjugate. The square 
brackets denote the separately gauge invariant pieces of 
each sector. As for the LO factorization, we use the Fierz 
identity to separate the spin and colour indices and match 
onto the appropriate jet and soft operators 



a (W6„,o) j(2 „6„) (rn) g j(o) {Tfi) g s (o ){Ts) 

+TC (la„,16„,0) j( o) (Tji) 3 j(0) (tr) S (0) (7>)< (43) 



The LO operators j(°>, and S^ ) are defined in Q 
and the subleading jet operator is 



J^) (/i;T) = 



Tr / da; + d d - 2 a;_Le 



^„ (0)iD£ (0) W,< 3 > (0, a£" ) (r)^ 3 ) (a^ , i fl )i£>x Q (z s ) Vn (**). 



(44) 




(a) 



T 

(b) 



(c) 



FIG. 3: Matrix element (0| J( 2q M |0) that give non-zero 
r values at 0(a s ). 



The operator is suppressed by A 2 due to the derivative in- 
sertions. The vacuum expectation value of this operator 
is shown in Fig. [3] 

In order to find the matching coefficients, we calcu- 
late the vacuum expectation value of (42) and (43). The 



Feynman rules for the dijet operators were written in [7] , 
and we find the vacuum expectation value of ( 42 ) is 



-a s r + 0(a 2 s ). 



(45) 



The vacuum expectation values of J^°\ and 
were found in ( 24 ) and the diagrams in Fig. [3] lead to 



(0\J^\fi,T)\0) = ±a s + O(a 2 s ). (46) 



Here and in the calculations below we have included the 
zero-bin procedure. As expected, the matrix element of 
the j( 2ob -) operator is suppressed by r ~ A 2 compared 
to the LO jet operator. The matching coefficients are 
found to be cf °"' lb "' 0) = -l + 0(a 3 ) and c£ an ' 1K ' 0) = 
+ 0(a 2 s ) meaning the hard functions are 



#l lQ "' lb "' V) = i + oK) 



H, 



(la„,16„,0) 



( A1 ) = + O(a 2 ). 



(47) 



Therefore, the contribution of (42) to the rate can be 



written in the factorized form ( 40 1 , with the appropriate 
subleading jet operator. 

B) i — 2d n , j — k — : The operator O^ 7 ^ is 
shown in ( |A4[ ) and accounts for the matching of QCD 
momentum conservation onto SCET momentum conser- 
vation at 0(A 2 ). The contribution where j = 2S n in 

(39) means taking O^ 6 "\o), which vanishes due to the 
explicit x dependence in the dijet operator. As in the 
previous example, we use the Fierz identity to factorize 
the contribution from this operator and match onto 



d d xOf K) \x)M {0) {T)O^ ] {Q) (48) 

= Cf 5 "'°'° 3 J^\r n ) ® J< m »>(t b ) ® S<°>(r.) 
+TCr A0) j(»>(r„) 8 j(°)(r s )®5(°)(r s ), 

where the new jet operators are 
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J(3«»)( T ) 



QN C 



Tr / dx + d d ~ 2 X±e 



l ^ + M0)W^H0,x^^M^(T)wi 3 \x^,x a ) (in ■ £ + in- d) (^)^(^) 



J( ° 5 " )(r) = ^ Tr / ^"^""^e-^"^-^^)^^^-)^^ (r)w| 3) (^% 0)^(0). 



(49) 



The derivatives in j( 25 ™) pull down the 0(X 2 ) component 
of the n-collinear momentum leading to a A 2 suppression 
compared to the LO jet operator. The x~ dependence 
of O^ 5 ^ is put into j( 015 ™) because only the n-collinear 
fields depend on this coordinate. We have introduced the 
Q's so the operators will have the correct mass dimen- 
sions. The diagrams of the vacuum expectation value of 
j( 2<5 n) anc i j(05 n ) are the same as the diagrams of the 
vacuum expectation value of and in Fig. [5J 

As before, the matching coefficient is most easily de- 
termined by c omp aring the vacuum expectation value of 



both sides of ( 48 1 . The explicit x in both O 



(25„) 



and 



j(o<5„)^ wn j c h is a derivative d/dp^ in mome ntu m space, 



acts on the (i-dependent prefactor in Ai^ of (|18j. There- 



fore, the vacuum expectation value of the left-hand side 



of (48) is 



a.r(---4Lj-3) +0{a 2 s ). 



(50) 



and the vacuum expectation value of the jet operators 
are 



^-2/,/-' r ) -Our) 



(0|j( 2 ^)(/.,r)|0) = a 
(0\J^(»,T)\0) = 25(T)+O(a s ) 



(51) 



Comparing both sides of ( 48 ) we find the matching coef- 
ficients are c{ 2S "' '°\^) = 1 + 0(a s ) and CQ 2S ""°'°\fi) — 
a s + 0(a 2 s ). We see in this example why the last line of 



( 39 1 is necessary: the d — 2 from the derivative acting 



on gives an extra term in (50) when it is analyti- 



cally continued to d = 4 compared to (0| J( os ™) |0), which 
has no poles at d = 4. The contribution of this operator 
product is thus factorized in the form of p0| ) with hard 
functions 



H 



(2<5„,0,0) 
1 

(25„ : 0,0) 



( M ) = -i + OK) 

in) = -d s + 0(a 2 s ). 



(52) 



Again, the matrix elements of the jet operators have the 
appropriate power counting and the logarithm in (51) is 
minimized at the jet scale, as expected. 

C) i — IS s , j — 0, k — In : As a final example, we 
show an insertion of the subleading measurement opera- 
tor. The contribution of the NLO measurement function, 
k = In, vanishes at 0(a s ) unless i — 1S S . This contri- 
bution accounts for the phase space we neglected at LO 
due to not conserving soft momentum. We factorize this 
operator product by matching 

' d d xOi Ua)f (x)M^ n HT)Oi 0) (0) 

= c[ 1S -°' ln} j(- 2 ^)(r n ) ® J^\rn) S^ m \t s ) 
+TC (Wn)j(0) {Tn) g j(o) (Tfi) S (0) (7i)< (53) 

The decoupling of the sectors in A^ 1 ™) means the 
NLO measurement operator can be treated identically 
as the LO measurement operator in (22) and gets pulled 



through when we use the Ficrz identity. Therefore, the 
appropriate jet and soft operators are 



./<-«•"»>(/,, r) 



I / dx+d d - 2 X ± e-^ x+ Xx^n^W^HO^X^^l^M^H^W^Hx^^n^niXn) 



N c 



Tr 



dt I dx-d d - 2 x ± e-^ x i> n (x n )W^>(x n) x™y^M%>(T)Wk a '(x%%nt)i> n (nt) (54) 



.(3) 



(0)/ 



r(3), 



QN, 



^&yW(a#\0) (iD Xa + S Xa ) ^Y^^O^^M^^Y^ix^^Y^^x^) 



where 



2fc (0)Q 

M^ a ( Tl {k}) = _±^x(o)( r ,{fc}). 



(55) 



becomes a derivative at infinit}Q We have also used the 



identity [3l |6] (in ■ d) 



= J dt i 



nt) to write 



the 1/pt, in (36) as a displacement in the j(° M1 ) opera- 
tor. The power counting suggests j(- 2S ^ M i) ^ \~ 2 w [\\ 



The explicit xj_ dependence in 0% " s/ is put into 
j{-25,M n ) b ecause we c hose the ft axis such that p R j_ = 0. 

The p*j_ in M„ n) acts as a total derivative at the cut SO 2 We use a co-variant gauge so the gauge field vanishes at infinity. 
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be enhanced compared to the LO jet operator. However, 
it is always convoluted with S^ s " Mn ' ^ meanm g the 
contribution to the rate will be 0(A 2 ) as expected. The 
diagrams for the vacuum expectation value of the jet and 
soft operators have the same picture as Fig. [2] and give 



(0|j(-^«)( M>T )| ) = 2d / (r) + O(a.) 
(0\J(°*"\»,T)\0)=8(T) + O{a s ) 
(0|S^ m »)( M ,t)|0) = -4a s r + 0(a 2 ), 



(56) 



which have the expected power counting in r. The 
matching coefficients are found to be c[ 1Ss '°' 1M "'\fi) = 



l+0(a s ) and C, 



(l<5 a ,0,lM„) 




(//) = 0+O(a 2 ) and the contri- 



bution from this operator can be written in the factorized 



form of ( 40 ) with hard functions 



H{ (ju) = 1 + 0[a s ) 



H, 







( M )=0 + O(a s 2 ). 



(57) 



This example shows how the decoupling of the measure- 
ment operators makes it straightforward to factorize their 
contribution. 

The factorization of the rest of the subleading dijet 
and measurement operators follows in the same way as 
the above examples. The explicit decoupling of the sec- 
tors makes finding the jet and soft operators a matter of 
using the Fierz identity to separate spinor and colour in- 
dices. The required jet and soft operators are written in 
Appendix [B] These operators are generalizations of the 
LO operators found in Sec. |III| The matching coefficients 
are found by comparing the vacuum expectation values 
of ( 39 ) . These matrix elements are pictured in Fig. [4] 
of Appendix [B] and their values are shown in Table [I] of 
the Appendix. The appropriate matching coefficients are 
shown in Table [TTJ. Combining the results of this Table, 
we find the 0(a s r) rate is 



R( T ) =l + a s ( -21n r-31nr+y -l + r(21nr-4) 



We can compare these results with those found using 
perturbative QCD in [13]. Although the results in [13] are 
for the exclusive two-jet rate using the JADE algorithm, 
the 0(a s ) QCD phase space is the same as the phase 
space for thrust, as mentioned above. Therefore, the full 
0(a s ) result calculated in [13] can be compared to (58 1. 



+0(a 2 s ,r 2 ). 



(58) 



Summing the rlnr's requires renormalizing the jet and 
soft operators, which we do not do here. 



V. CONCLUSION 

We have shown how to systematically calculate the 
0(t) corrections to the thrust rate in the perturbative 
regime using SCET. The rate was factorized and writ- 
ten as the convolution of the vacuum expectation value 
of jet and soft operators. Each operator has consistent 
power counting and depends on a different scale asso- 
ciated with the rate. The appropriate jet and soft op- 
erators are found in this work, as well as the matching 
coefficients. The 0(a s r) rate was calculated and repro- 
duced the rate found using perturbative QCD. 

The rate was factorized by matching the QCD cur- 
rents and measurement operators onto SCET dijet and 
measurement operators, which in the formulation of [7J, 
explicitly decouple the n-collinear, n-collinear, and soft 
degrees of freedom. The non-local product of the dijet 
and measurement operators were then matched onto jet 
and soft operators that separately describe each of these 
degrees of freedom. The approach illustrated here can 
be applied to other jet observables to calculate sublead- 
ing corrections. We are currently exploring subleading 
corrections to the continuous angularity observables of 
which thrust is an example. 
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Appendix A: Dijet Operators 

The LO dijet operator was given in Q. In this section we write the NLO and N 2 LO dijet operators necessary for 
calculating the 0(a s r) thrust rate. They are found by following the SCET formulation of [7J. The NLO corrections 
to the n-collinear sector are described by the operators [7J 



O 



(x) 



0^\x) 



^{xn^Up^x^TW^Xx^xf) 



( X ™° , )yi 3) (0, x™ " ) Wi 3) (x™ ,x n )i> n (x n ) 



,(3) 



^O^n %'-0^'n (Oj X n ) W n (x n , X n )P n 1p n ( Xn ) 
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Oi le "'{x) =-i I dt 
Jo 



^ n (s« + tn)T d ^ lvl TW^{xn + tn, x Tl )P n ij n {xn)W^ dc {x Tl + tn,x°°) 



(Al) 



o£ M {x) = -i I dt 



^n{x n )P n T lvl ^W^\xn,x n + tn)T d ^ n (xn + tn)W^ dc (x n + tn, x% 
Y^x™* , 0)Fi 8) >, x™')] [n M wj? ) *{x%,x n )igG?r{x n j 

The covariant derivatives are defined as = — igA^ and the gauge fields are in the same sector as the fields the 
derivative is acting on. The m th -sector field strength tensor denoted as Gif% . For the corrections to the n-collinear 
sector, the two operators required are [7] 



0^\x)= \Mx n )wi a Hxn,xT)} \YW(x™°,0)Yi 3) (0,x°°') 



W^ 8) (x~,s„)ri^x(a; n )|v ft (a n ) 



(A2) 



See [7] for the Feynman rules and physical pictures of these operators. The matching coefficients are 

=cf as) = -l + 0(a 8 ), Ci lbn) =cf M =l + 0(a s ), cf Bn) = -1 + 0(a s ), cf M = 1 + 0(a«). 

(A3) 

The operators describing the N 2 LO corrections to the n-collinear sector are found in the same way as the NLO 
operators. The operators are 

0f an \x)= \^ n {x n )P n Tin- D(x n )Wi 3 \x n ,x^)] [yW(ar~',Q)r B (a) (0,^)1 M 8) (a:~, i n )^4(in) 



{ 2 bn) (x) = -i I dt 



i) n (x n + tn)$) ± (x n + tn)^T^Wi 3 \x n + tn, x n )i]p ± (x n )W^ (x n , x%) 
Y& , 0)yi 3) (0, a£" )1 [^ 3) {x™,xMx n ) 



Of An \x) — i / ^ n (i fi )P fi wW(a: fi ,a:fl + tn^^Gf + tn)W< a) (a; fi + tn,x?) 



Y^ (x%' , 0)ri 3) (0, x°°° )} \w (3) (x%,x n )P n M*nj 



(A4) 



o 



(2S r 



\x) = Q ^ n {x n )P n T (n-D + n- (^''(is,^)] [yW(x™°,0)Y± 3 \o,x%>°)\ [x~wf' \x% ,x n )P^ n (x n 



with matching coefficients 



a 



f an) =-\+0(a s ), cf K) = -1 + 0(a s ), C? An) = l + 0(a s ), C^ = l + 0(a s ). (A5) 



The Q in 0^ Sn ^ is required because of the definition in The operators describing corrections to the n-collinear 
sector are 

0< 2aft) (x) = [^(x^PnW^Hx^x^)] [Y^{x^',0)Y^\Q,x^-)] \wi s \x^,x n )Tin-D{x n )P n ^ n {x n ) 
0^\x) = -i dt Mx ri )Wi 3 \x n ,x^) yW(x™%0)y( 3) (0,x°°° 



Wi 3 \x™,x n )t1p ± {x n )wi 3 \x n ,x n + tn)^Tpp(x n + tn)^ n {x n + tn) 



cf =-i I dt 



jn(xn)PnWW(xn, x°?)\ [F„ (3) «% 0)y fi (3) (0, xT) 
W^\x™,x n +tn)gGf(x n + tn)<j ±a0 rP n Wi 3 \x n +tn, x n )M*n) 



(A6) 



0<£ &a) {x) = Q [x + ^ n (x n )P n w( 3 Hxn,x°°)\ [ri 3 )(^%0)yi 3) (0,x~ s )] [wi 3 \x™,x n )T (n-D + fi-ty (x n )P n ip n (x n ) 
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with matching coefficients the same as their n-collinear counterparts. The intuitive picture of these operators is similar 
to the pictures presented in [7]. 

The NLO corrections to the soft sector are described by the operators [7] 



0^\x) 



i I dt 
lo 



O 



(lc fio ) 



Y^\x^,tn)i5 ±fi (tn)Y^(tn,0)YP(0,x^)\ [wf \x^,x n )PM^) 

/>00 

$ n { Xfl )PnTWW (a*, x%)\ (x™° , 0)Y± 3) (0, tn)iD Xll {tn)Y^ ] (tn, x°?' ) 
w£\x™,x n )iDl_(x n )Pnipn(xn) 



(x) = —i I dt 
Jo 



(A7) 



0$*'\x) = Q \x^ n {x n )P n YW£\x n ,x™)\ \Y( s \x™>,0)(D» J _+fy x )Y« s \0,x?°°)\ Iw^ ] (x™ ,x n )P n ^ n (x n ) 



0%*»\x) = -i 



dt 



n^gGT{x n )W^ ab {x^xt) 

i 



r„ (8)fcc (x^, tn)A (tn)T%^ (tn, 0) '^ ± TY^ (0, x°°° ) 



i I dt 
lo 



MXn)PnWi 3) (Xn,X°?) 

wi 8)ba (x™,x n )n^gGT(x n ) 



Wi 3 \x™,X n )PMXn) 



Y^ , o)r 7l/± |r c ri 3) (o, tn)^ s (tn)Y^ )ch (tn, x% 



with matching coefficients 

= C { 2 lcris) = -2 + 0(a s ), C^ ldns) = Of = 1 + 0{a s ), C [ 2 15s) = 1 + 0(a s ). 
The Feynman rules of these operators are given in [7]. The N 2 LO corrections for soft gluon emission are 



(A8) 



O. 



(2A„„) 



(x) = -i 



dt 



Mx n )P n Wi 3 \x n ,x°°)} \Y^(x^,tn)ga^G^(tn)Y^(tn,0)TY^ 3) (0,x^) 



{ 2 Ms) (x) = -i I <ll 



X Wl 3) (x™,X n )P n i> n (Xn) 

with matching coefficients 



r (2A n .) _ ^(2A flfl ) 



1 



+ 0(a,). 



y2 _ " 2 2 

The corrections from the expansion of the momentum conserving delta function are described by 



(A9) 



(A10) 



of- 


(x) = Q' 




{x) = Q 


Q (2S S± 


w-f 




X 



(ar) = Q [^„(a; fi )P fl rW ) l 3) (^,^)] [V^ (a:~' , 0) (n- P + n- y B (a) (0,a£")] [wi 3) (x™,x n )P n ip n (x n ) 
(x) = Q ^ n (xn)P n TWi 3 Hxn,x°°)\ [^)(C,0) (n ■ D + n ■ ri 3) (0,a;^ s )] [aT wf } (a£\ x„)P fi Vn(^) 



^^^(x^p^r^i 3 )^,^)! [ri 3 )«\o) (p^ + ^i^i) r« (3) (o,a£ 



W n (^yj, 7 ^n)-^nVVi (^n) 



o. 



(2<5c 3 



'(a;) = -i / di 



^^(^jPfiK^K^)^ 3 ^^,^) 



(All) 



ri 3 )(^,M^^(^)^ (3) (^,0)(P^ + ^)(0)Fi 3) (0,^)1 fw^ 8) (x~,a n )P fi ^ a (a; n ) 



0^(1) = -i 



^(a^P^ 3 ^^) 



x yi"(C*,0)K + P^)(0)ri 3) (0,^)iP^(M)ri 3) (tn,a;^) wf ) (^^„)iP l/± (a: n )rP fl 7/; fl (x„) 
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(a) (0|J<°8)|0> (b) <0|J( 2£ )|0) (c) (0|St 2d »»)]0) (d) (0|S(°8)|0> 



FIG. 4: Diagrams for the operators describing a soft quark and the quark and anti-quark in the same sector. The 
double dashed lines represent Wilson lines in the adjoint representation. Reverse the fermion arrows for (0| J^ 2 ^'|0). 
The filled boxes highlight a Wilson line changing representation. 



with matching coefficients 

Cf 5 ° +) = C< 2 tf -> = 1 + 0(a s ), c¥ Ss±) = H-0(a»), C 2 (Mc * n) =C^ } = -2 + 0(a s ). (A12) 

We do not require the N 2 LO corrections from soft quark emission for an N 2 LO calculation. 

As expected, the subleading operators can be written as separate n-collinear, n-collinear, and soft pieces. The 
Feynman rules for the LO and NLO operators at 0(a s ) were shown in [7]. For the sake of brevity, we do not show 
the N 2 LO Feynman rules here. We are only concerned with vector currents in the above calculations, so T = 7^. 



Appendix B: Jet and Soft Operators 



The jet and soft operators are found by doing the matching in (|39|. It is convenient to use the basis 



1 Q ctB 



2 ' 2 



i ft i ft i * 
1 7j_75? 2' 2' 2 2 2 2 



(Bl) 



for the Dirac matrices. We have defined 7^ = g^j v and = | [7^ , 7^] where 5^ = — (n^n v + n"n^)/2. For use 

,Jh v t a ^ >lv . The jet and soft operators will be parity even scalars due to only considering 



a/3 



later, we aslo define e 
vector currents. 

The leading order jet operators required at 0(a s ) are 



J (08 W) 

j (05a W) 

j(")( M ,r) 



1 

5j_ 7 



Tr / dx + d d - 2 x ± 



QN C 
Q_ 
N c 
1 



N c 



dx+d d - 2 x ±e - lQx Gr^)^ (O-C^rW^' to (^, Ifl )Gf(s R ) 
Tr / ^+d d -V e -^V^(0)W^^ 



Tr / dx + d d - 2 x ± e 



Although J(° c "= 5 ») has a £> ± suppression, there is an enhancement by an explicit x± so the operator is LO. The 
vacuum expectation value of the operator j(° 8 ) is shown in Fig. 4a The other operators have the same diagrams as 
Fig. [2] with the appropriate derivative insertions. The Q's are introduced by dimensional analysis. There are also the 
enhanced jet operators 



/- m %t) = |^Tr / ,/,--,/'' 2 r_< 



■«* + (x i ) J i(0)^(0,C)k 0) W^ ) (^^)^W (B2) 



N c 



2 Q dr" 



j(- 2M ^)( M ,r) = — Tr / dx+d d - 2 i r ±e - iC 3 a:+ ^ n (0)^ 3 )(0, a; ^) 



fin ■ d d 
2 Q ~£h 



M^{r)W^\x^,x n )^ n {x n ) 
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I 


<O|J«0*,t)|O) 





5(t) 


8 


-8(t) 


0<5„ 


-2S(t) 


Oc n3 5 3 


-25{t) 


-25 s 


25' (t) 


-2S s M n 


S'(t) 


-2M nb 


S'(t) 


2ab n 


d s /2 


2a n 


a s (-2/e-2Lj-9/2) 


2A n 


d a 


2S n 


a s (-2/e-2Lj-3/2) 


2e 


d s (-2/e- 2Lj + 2) 


2/ 


a s (-2/e-2ij+2) 










5(t) 


8 


-J(t) 


Omi 


S(r) 


UM2 


d(r) 


0<5„ 


-2S(t) 


-2M Sl 


S'(t) 


2afefj 


ds/2 


2a^ 


a s (-2/e-2Lj-9/2) 


2A.fi 


a s 


2Sn 


a s (-2/e-2Lj-3/2) 



(b) 



(a) 



n 


/nl <?(")('// -r"! ln\ 





5(t) 


8 


5(t) 


2c ns $s 


a s (l/e + L s - 1) 


28 C ns 


a s (l/e + L s - 1) 


2dns 


a s (-l/e-L s + 1) 


2dfis 


a 8 (-l/e-L s + 1) 


2A ns 


d s (2/e + 2L s -2) 


28 sn 


a s (2/e + 2Ls -2) 


2M S1 


-4a s 


2M a+t 


— 4« s 


4<S a x 


— 4a s r 


45 3 M„ 


4a s r 


4M„ 6 


2a s r 


4M nb 


-2d s r 



(c) 



TABLE I: Relevant vacuum expectation values of the jet and soft operators for the 0(a s r) thrust rate, 
operators distinguished by primes give the same values. 



The 



where the derivatives act at the cut. These operators are always associated with 0(A 4 ) soft operators so do not give 
an enhancement to the thrust rate. We have neglected the contribution from .A/f(~ 2 ™") because it vanishes at 0(a s ). 
The required N 2 LO jet operators are 

J (2aK) {^r) = ^-Tr J dx+d d - 2 x±e-^ + AMDl{mk 3 \^Xn°)*M^ 

J (2e W) = -^Tr f X ' dsdt [ dx+^ixe^Vf'^x^^ + ttlitx^fti^. + tt) 
Q^c Jo J 

J (2/) (^,t) = -^—Tr / dsdt [ dx + d d - 2 x±e- lQx+ Wi s)a \x™,Xn + tn)4> n (xn+tn)T b W^ 
Q^c Jo J 

x \ ll^Mxn)M£^ n (0) W±a ^Wi 3 \0, sn)T l ^ n {sn)W^ la {tn, x™*) (B3) 
J (2a "W) = g^Tr| dx+d d - 2 x^e-^ + MQ)in- D{Q)W^\Q 1 x^)^M^ 
J (2A -H^t) = q^KJ dt J dx+d d - 2 x ± e-^ + M0)Wi 3 \0M9^Gn P (tmi 3 Htn,,x^) 

^l,M^{T)W^{x™,x n )Mxn) 
J (25 "W) = ^-TtJ dx + d d - 2 x^e-^ + MO)Wi 3 \0,x^)^M^\T)W^ (in t + in ■ D) (x n )M^n) 

and their hermitian conjugates. The vacuum expectation value of j( 2 ° b ") was pictured in Fig. [3] Although J< 2e ' and 
j( 2 ^ look complicated, they have simple diagrams shown in Fig. 4b The 75 in j( 2j4 ™) is necessary because e™' 9 is 
parity odd. The vacuum expectation value of these operators are shown in Table [la| 
The J operators are similar to the J operators. The leading order operators are 

J ( °W) = j^-Tr J dx-d d - 2 x^e-^^{x n )w!i i \x n ,x™) 1 tMf{T)W^ 

j(Os) (/i)T) = afn^u I ^-^J^-^^^^^^x^^^Jia^^j^a^ 
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jCOmO^) = ^Tr^ dt J dx~d d - 2 x±e-^ Mx n )wl?\x n ,x™)^Mf {t)W { ^ ] {x^ ,tn)^ n {tn) 

J (05 "W) = ^Tr /" d,-rf d - 2 Ii e- ! « I ^^ s (.„)^ 3) ( I „,C)^^r)W'i 3) (x^,O)^(0). 
There is only one enhanced operator that needs to be considered 

j(- 2 ^)( M ,r) = ^-Tr| dx-rf'-VLe-^fiW^ 

(B4) 

The N 2 LO operators are 

J^)(m,t) = — ^Tr / dx-d d - 2 x±e- i Q*-$ n {x n )w£\x n ,x™»)t ) M%\T)w£\x™ (B5) 

j( 2 ^V,r) = JL-Tr^ 00 ^ | &+d d - 2 xxe- iC2a;+ Vi„(x„)^ 3) ( a;n ,07 5 |Alf (r) 

x Wi 3) (x°° , tn)gef Gf (tn) (tn, 0)^ (0) 
J {2S *Hfi,T) = ^-Tr J dx-d d - 2 x ± e-^Mx n ) (in-*D+in-D) (x n )w£\x n , x™>)^M.f (t)w£\x™ , 0)^(0) 

The diagrams are similar to those found for the J operators. For example, j(° 8 ^ is the horizontal reflection of Fig 



4a 



lb The difference in the required J operators 



The vacuum expectation values of these operators are shown in Table 
compared to the J operators is because we have chosen the n axis to be anti-parallel to the n-collincar sector. 
The leading order soft operators are 

S (0) (M, t) = -^Tr F fi (3) (x°° ■ , 0)Fi 3 > (0, x™* (r)Y^ (x™* , 0)Y^ (0, x?- ) 

iVc 



S ( ° 8 W) = ^r fi (8)o ^%0)yW 6c (0,^ (B6) 
which only differ in the representation of the Wilson lines. The vacuum expectation value of S^ 08 ' is pictured in Fig. 



4d The N LO operators are 



S(2c„A) (M)T ) = _L_ Tt J dt y fi (3) (a; oo %0) (fia + ( 0)Y ( 3 \0,x^)M^(t) 

x yj 3 > (x™° , *n)i^5 (tn)yi 3 ) (tn, 0)Y^ 3) (0, a#* ) 
S^\u,t) = —TtJ^dtY^^^Q) (Sl + iDl) (0)Y( 3 \0,tn)D ±a (tn)Y^(tn,x^)M^(r) 

S^\a lT ) = dsdtY^ a \x^,sn)Usn)T b Y^ 

xY^ 3 \x^%0)Yj l 3 H0,tn)T c Mtn)Y r i 8 '> ca (tn,x^) 
S^(a,r) = -^Tr ^ 'dsdtyi^^sn)^^ 

xri 3 )(^ % o)yi 3 Ho,in)r^sMri 8)ca (in,^=) (B7) 

— ^oo 

5 (2A„ s)(M)T) = 9^l Tr / A1 K») (l «. )0)y W (0)to) ^ (te) y( 8 ) ( ^ C) ^0) (t) ^3) (C 
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Q-^c Jo 



1 



TrK 



(3) 



QN C 



(a£" , 0) (in • £ + in • D) (O)Y^ (0 , x %. )Aff (r)Y^ (x%- , 0)F fi (3) (0, < ' ) 
^. .. „ (<\0) («•£ + «•£>) (QJ^HO.^O^W^W'-O)^^^') 



There are also N 2 LO operators involving the measurement function expansion that we write as 
S (2Mi) (m, r) = ^Tr yf } (a-- , 0)F„( 3 ) (0, )A?f > (t)Y™ «>* , 0)Y ^ (0, x% 



(B8) 



where 



2f? (0) 'i 2 f> 
M^(r,{k}) = P -^^^M^(r,{k}) 



(B9) 



where the cancels in the definition of k\ t in ( 32 ) . The contribution from the other measurement operators in ( 37 ) 
vanish at 0(a s ). 

There are also N 4 LO operators that contribute to the N 2 LO rate due to convoluting with enhanced jet operators. 
These operators are 

S^(^r) = J-Try fi (3) (:r-\0) fi&±a + DJ_D ±a ) Y™(0, x^)M^ (r)Y^(x^, 0)Y^ (0, x%°) 

far) = -±-TrYW(xZ>-,0)(iDZ + l ^l)Y^ 3 \o,x^)M^ a (r)Yj i 3 \x^ ,0)Y^(0,x^) (BIO) 
S^\fi, r) = -^Tr ri 3) (*?>* , 0)yj 3 ) (0, x™* (t)Y^ ( x %. , )y fi (3) (0, ^ s ) 



where 



-ffc (0) ) 2 

M^(r,{k}) = ^ ^< 0) (r,{fc}) 



(k {0) ) 2 

%^-Mi 0) (r,{fc}). 



(Bll) 



The vacuum expectation value of the soft operators are shown in Table [TcJ 

The operators in Appendix [X] are matched onto combinations of the operators in this section. The matching 
coefficients are found by taking the vacuum expectation value of both sides of ( 39 1 and are shown in Table [TTJ Using 
this Table and Table [I] the 0(a s r) rate in (58) can be calculated. 
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(i,3, k) 


(I, Tn, n) 






R(r) 


(lfln, 16n, 0) 


(2a6„,0,0) 


-1 


+ 1 


a s r/2 


lo„,0) 


(2afct ,0,0) 


-1 


+ 1 


a s r/2 


(2o n ,0,0) 


(2a„,0,0) 


1 


-1/2 


a s r(-l/e-Lj-9/4) 


(0,2a n ,0) 


(24,0,0) 


1 


-1/2 


a s -r(-l/e - Lj - 9/4) 


(2A„,0,0) 


(2A ni 0,0) 


-1/2 


-1/4 


-a s r/4 


(0,2A„,0) 


(2Aj,,0,0) 


-1/2 


-1/4 


-a s r/4 


(2<5 n ,0,0) 


(2<5„,0<5„,0) 


1 


1 


& s (4/e + ALj + 7) 


(le n , le„,0) 


(2e,0 8 ,0 8 ) 


-1/2 


-1/2 


a s r(-l/e- Lj -2) 


(l/ n ,l/ n ,0) 


(2/, 8 ,0 8 


-1/2 


-1/2 


a 3 T(-l/e- Lj -2) 


(lo ft) 16 ft) 0) 


(0,2a6 n ,0) 


-1 


+ 1 


a s r/2 


(lferi, lafi, 0) 


(0,2o6t,0) 


-1 


+ 1 


a s r/2 


(2o ft ,0,0) 


(0,2a ft ,0) 


1 


-1/2 


a s r(-l/e- Lj -9/4) 


(0,2a n ,0) 


(0,24,0) 


1 


-1/2 


a s r(-l/e-Lj-9/4) 


(2A B ,0,0) 


(0,2^,0) 


-1/2 


-1/4 


-a s r/4 


(0,2A n ,0) 


(0,2,4,0) 


-1/2 


-1/4 


-a a r/4 


(2<5 n ,0,0) 


(25„,0<5„,0) 


1 


1 


a s (4/e + 4Lj + 7) 


(l<5 s ,lc ns ,0) 


(OCnsSs, 0, 2c ns Ss) 


1/2 


1 


a s r(-2/e-2L s -2) 


(2<5c sn ,0,0) 


(0c„ s S ai 0, 2<5c ns ) 


1/2 


1 


a s r(-2/e-2L s -2) 


(ldns, ld„ a ,0) 


8 ,0,2d ns ) 


1 


1 


a s r(l/e + L s + 1) 


(ld fts , ld ns ,0) 


(0 8 ,0,2d Ss ) 


1 


1 


a s r(l/e + L s + 1) 


(2A„ S ,0,0) 


(0,0,2A„ S ) 


1 


1/2 


a s T(l/e + I s + 1) 


(0,2A ns ,0) 


(0,0, 2Al s ) 


1 


1/2 


a s r(l/e + L s + 1) 


(2A„ S ,0,0) 


(0,0,24„ s ) 


1 


1/2 


a s r(l/e + L s + 1) 


(0,2^,0) 


(0,0, 2AL) 


1 


1/2 


a s r(l/e + L s + 1) 


(25 s _,0,0) 


(0,05„,25 s _) 


1 


1 


a s r(-4/e - 4L S - 4) 


(25 s± ,0,0) 


(-2<5 S ,0,4<5 S± ) 


1/2 


1/2 


4a s r 


(15„,0,ln) 


(-2«y a M n ,0 M i,4J 8 M„) 


1 


1 


4a a r 


(0,0, 2n b ) 


(-2M 2nt ,,0M2,4M„ b ) 


1 


1 


— 2a s r 


(0,0, 2n 6 ) 


(0,-2M ftb ,4M ni> ) 


1 


1 


2a s r 


(0,0,2*0 


(0,0 M ,2M s i) 


1 


1 


-4a s r 


(0,0, 2*+*) 


(0,0m, 2M +t ) 


1 


1 


— 4a s r 



(j,J,fe) 








(2<5„,0,0) 


-4a s 


-4a s 


—4a 3 T 


(le„, le„, 0) 


a s 


a s 


a s T 


(l/n,l/n,0) 


a s 


a s 


a s r 


(2^,0,0) 


-4a s 


-4a s 


—4a s T 


(ld ns , ld„ s , 0) 


— 2 S 




~a s r 


(ldn S , ldns, 0) 


— 


-Qs 


—a s T 


(2<5 S _,0,0) 


4a s 


4a s 


4a s r 



TABLE II: The operators and matching coefficients for the 0(a s ) factorization. The dagger means the operator is 
the hermitian conjugate. At the order we are concerned with, there is at most one set of subleading jet and soft 
operators to be matched onto, so a = 1. The table at the bottom gives the matching coefficients for the last line of 



(39 1. Values of the hard functions, which are defined as H [ ^' k> = C^C^'C^ q' " , are also given, as well as the 



contribution of these operators to R(t). 
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